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Abstract. The orbital eccentricity of a single planet in a stellar binary sys- 
tem with a sufficiently large mutual inclination angle is known to oscillate on a 
secular timescale through the Kozai mechanism. We have investigated the effects 
of the Kozai mechanism on double-planet systems in binaries. The evolutionary 
sequence of a pair of planets under the influence of a binary companion is fairly 
complex. Various dynamical outcomes are seen in numerical simulations. One 
interesting outcome is the rigid rotation of the planetary orbits in which the 
planetary orbital planes secularly precess in concert, while the orbital eccentric- 
ities oscillate synchronously. In such cases the outer planet acts as a propagator 
of the perturbation from the binary companion to the inner planet and drives 
the inner planetary orbit to precess at a rate faster than what is predicted by 
the Kozai mechanism. 



1. The Kozai Mechanism 

The first decade of radial- velocity and photometric surveys has revealed that 
planetary systems commonly form among stellar binary systems. As of July 
2007, at least 37 planetary system s (~ 20% of all the known systems) are 
mem bers of multiple-stellar systems ( Raghavan et al.ll2006l : iDesidera fc Barbieril 



2007). The true stellar multiplicity among planetary systems may be higher 



than 20% since the photometric surveys in search for stellar companions around 
currently known planetary systems are not yet complete. Because most of these 
stellar companions are in wide orbits (typically a = 10 2 AU or greater), their 
orbital inclinations are not correlated with the invariable plane of the planetary 
systems, implying that many binary companions orbit around planetary systems 
in largely inclined orbits Q. 

It has been known in celestial mechanics that in a hierarchical triple-body 
system with sufficiently large mutual inclination angle (I > J cr it = 39.2°) be- 
tween the inner and outer binaries, the inner binary's orbital eccentricity and 
inclination undergo larg e- amplitude oscillations on a secular timescale through 
the "Kozai mechanism" ( Kozai|[l962l ). When adopted to planetary systems, the 



inner binary orbit can be replaced by the osculating orbit of a planet around the 
primary star. The Kozai mechanism has following unique properties (see also 
Figure [T]): (1) the orbital eccentricity of the planet oscillates with an amplitude 



If the distribution of the mutual inclination angle is completely isotropic, the probability of 
the mutual inclination angle being greater than the critical Kozai angle (39.2°) is 77% 
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Figure 1. An example of secular evolution of a planet orbit through the 
Kozai mechanism. The amplitude of the eccentricity oscillation (e max ) is de- 
termined by the initial relative inclination angle. The argument of periastron 
typically alternates between libration and circulation during the Kozai cycles. 



that is dependent (in the lowest order) only on the initial mutual inclination 



angle Iq; e max = y 1 — 5/3 cos Iq ■ Other orbital parameters only affect the 
eccentricity oscillation period, 

Pe^/PplU-^) 3 / 2 , (1) 



( Ford. Kozinskv <fe Rasid 120001 ) where P p \ and are the orbital periods of 



the planet and the binary, respectively, and eb is the binary eccentricity. (2) 
The orbital inclination / of the planet with respect to the binary plane also 
oscillates on the same timescale as the eccentricity oscillation, conserving the 
Kozai integral Hk = \/l — e 2 cos /. (3) The argument of the periastron u al- 
ternates between the librat ion mode around ±90° and the circulation mode 
dFabrvckv & Tremaindl2007l ). 



An episode of the Kozai mechanism in an extrasolar system may leave a 
directly observable signature in the orbital eccentricity of the planet. Theo- 
retical models show that several highest p lanetary eccentricities observed are 
best explained by t he Kozai mechanism ( Holman. Touma &: Tremaind 1997 ; 



Wu &: Murray 20031 ). Also, the median eccentricity of the planets in bina- 



ries (eme^binary = 0.25) is higher than that of the planets around single stars 
( e med,singie = 0.20), implying that the planetary eccentricities are commonly 
excited by stellar companions. 

The Kozai eccentricity oscillation results from the torque exerted by the 
binary companion averaged over the planetary orbit. If there is another secular 
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perturbation at play that precesses the planetary orbit in a shorter timescale 
than the Kozai oscillation timescale (Eq. [1]), the averaged Kozai torque is re- 
duced, and the Kozai eccentricity oscillation will be suppressed. Extra preces- 
sions are typically caused by general relativistic effects, tidal torques, or m utual 
interactions between planets ( Wu Murray! 120031 ; iTakeda Ras"ioll2005l ). 



2. Evolution of Double-planet Systems in Binaries 

2.1. The Classical Secular Theory Multiple-Planet Systems 

Although the Kozai mechanism has been successful in modeling single planets 
in binaries, there are only handful of literatures that investigated the evolution 
of multiple-planet systems in binaries. What makes the problem challenging is 
the large dimensions of the parameter space; six orbital elements and the mass 
for the two planets and the binary companion already sum up to 21 parame- 
ters, virtually prohibiting a simple numerical exploration of the entire parameter 
space. Analytical understanding is necessary to identify different evolutionary 
outcomes. 

The extra complexity for the evolution of double-planet systems in binaries 
is due to the additional secular torque which planets exert on each other in addi- 
tion to the Kozai torque applied by the binary companion. All the torques result 
in the secular precession of planetary orbits and evolution of orbital eccentrici- 
ties. The secular orbital solutions of planetary systems through mutual planet- 
planet interactions have been formulated through the classical Laplace-Lagrange 
theory (L-L theory, hereafter). Assuming small initial orbital eccentricities and 
inclinations 0, the L-L theory can predict from initial conditions the evolutionary 
sequence of the orbital parameters, namely e(t),Lo(t), I(t), and Q(i). It should 
be kept in mind that the L-L theory is an approximated perturbation theory 
and is not a reliable estimator for predicting accurate long-term behavior of a 
planetary system. However, the great merit of the L-L theory is that one can 
extract useful evolutionary information through simple linear algebra and derive 
various order-of-magnitude estimates for evolutionary timescales. 

2.2. Numerical Simulations 

We have run a large set of numerical simulations with various initial conditions 
to investigate different evolutionary outcomes of hierarchical four-body systems 
(i.e., double planets orbiting a component of a wide stellar binary). The equa- 
tions of motion were fully integrated by the Bulirsch-Stoer method to account for 
all the possible close encounters between planets or between a planet and the pri- 
mary star. To reduce the number of free parameters, we have used fixed binary 
parameters: a-Q = 750 AU, Ib = 50° (with respect to the invariable plane of the 
planets), ee = 0.2, and both binary components with masses Mq = Mb = 1M Q . 
Also, initially the two planets are both in near-circular orbits with eccentricities 
less than 0.01 and separated by more than one mutual Hill radius to ensure 



2 Typically, the terms up to the second order in e and / are included in the st andard L-L theory. 
The g eneralized fourth-order L-L theory has been recently developed by IVeras fc Armitagel 
(|2007h . 
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Figure 2. A numerical simulated evolution of rigidly rotating double-planet 
system. Dashed curves represent the orbital elements of the outer planet. All 
the angular quantities are denoted in degrees. The longitudes of ascending 
node and the orbital inclinations of the two planets are almost indistinguish- 
able throughout the evolution. The orbital eccentricities of both planets oscil- 
late at the same period. Notice that when e-i is excited to ;> 0.1 by the Kozai 
mechanism, the periastron arguments get locked into tight apsidal libration. 



that the system does not immediately become unstable. All the systems were 
evolved for 500 Myr, a time sufficient for the outer planet to undergo at least a 
few Kozai cycles, unless the system has gone unstable. 

2.3. Rigid Rotation of Planetary Systems 

A variety of dynamical outcomes have been observed in our simulations. Of 
particular interest was the rigid rotation of the planetary system, in which the 
orbital elements of both planets secularly evol v e in concert. A similar result 
has been previously observed by llnnanen et all ( 19971 ). They have numerically 



calculated the evolution of the Solar planets with the presence of a hypothetical 
stellar companion to the Sun and observed that the inclinations and ascending 
nodes of the planets are always aligned throughout the evolution, as if the planets 
had been embedded in a rigid disk. 

Figure [2] shows an example of a rigidly rotating double-planet system. It 
is immediately noticeable that the orbital inclinations and the longitudes of 
ascending node are strongly aligned throughout the evolution, maintaining me- 
chanical rigidity of the orbital planes. Another interesting aspect of Figure [2] is 
that ei oscillates at the same periodicity as but with a smaller amplitude. 
These synchronous eccentricity oscillations are contrary to what would be pre- 
dicted by the Kozai mechanism; a pair of coplanar planets at different orbital 
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Figure 3. Secular contours drawn in the Auj - e2/ei space, from a numerical 
simulation (left) and the classical L-L theory (right). The numerical simula- 
tion started from initially equal eccentricities of 0.01, with mi = 0.05Mj at 
1 AU and — lMj at 8AU. Each cross is 1,000 years apart, starting from 
the circulation regime at e^je-i — 1 and Au>o = 90. In ~ 0.3 Myr, the system 
evolves past the secular separatrix (dashed curve) and eventually falls in to 
the small amplitude libration island in ~ 10 Myr. Similar trajectories derived 
from the L-L theory are drawn in the right panel for comparison. 



radii are expected to undergo eccentricity oscillations at equal amplitudes but 
with different periods. 

The synchronized eccentricity oscillations result from the small-amplitude 
apsidal libration, for which Aui = \u>\ — o^J < $ is maintained fo r a small angle 
<F It has been previously discovered bv lChiang &: Murray ( 20021 ) that the adia- 
batic increase of the outer planet's eccentricity due to the remnant gas disk drives 
the apses of the planets from independent precessions (circulation) to strongly 
coupled precessions (aligned libration) . The same apsidal trapping mechanism is 
essentially what causes the Kozai-induced synchronized eccentricity oscillations. 

The numerical result of the sample system in Figure [2] is plotted in an 62/ e± 
- Aw phase diagram in Figure [3] (a). Figure [3] (b) shows the secular contours 
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computed from the L-L theory for the same system, with varying initial eccen- 
tricities and periastron arguments. The secular separatrix in the Figure [3] (a) 
and (b) is drawn from the apsidal libration criteria such that 

e 2 n 5 qa^ [l - (l/8)a 2 ] 

< ~o 1 T72 cosAwo (2) 

e± t o 2 1 — qa 1 ' z 

and 

2 1 - qa 1 / 2 1 
ei, 5a[l - (l/8)a 2 ] cosA^o 



( Laughlin. Chambers h Fischerll2002l : IZhou Sunll2003l ). Here, e\ : o and e2,o are 



the initial eccentricities of the inner and outer planets, respectively, q = mi/m,2, 
a = ai/02, and Aloq = uo^ — cui t n is the initi al angular offset between the apsidal 
vectors of the planets. IZhou Sunl ( 20031 ) called the libration modes defined 



by Equation ([2]) and ([3]) "down-libration" and "up-libration" , respectively. The 
separatrix in Figure [3] corresponds to the up-libration boundary for which the 
two apsidal vectors are aligned (they are anti-aligned for the down-libration). 

The evolutionary paths of two planets around a single star can be deter- 
mined solely by their initial conditions, i.e., for each set of initial orbital param- 
eters there is one corresponding secular contour that dictates the evolutionary 
sequence of the planets. With the presence of a binary companion, however, the 
outer planet's eccentricity is secularly excited through the Kozai mechanism, 
permitting the system to explore different evolutionary modes 0- 

In Figure [3] (a), the simulation starts at eo,i/eo,2 = 1 and Aujo = 60°. The 
system is initially placed in the circulation regime so that the two apses precess 
independently. However, as e2 starts to evolve significantly through the Kozai 
mechanism while e% remaining at a small value, the eccentricity ratio increases, 
driving the system rightward in the e2/ei - Auj space. The system marches on 
eventually past the secular separatrix, entering the libration regime. Further 
increase in e2 brings the system toward progressively smaller libration islands 
and therefore damps the libration amplitude. The system maintains a small- 
amplitude secular libration while e2 is sufficiently large (;> 0.1), until the outer 
planet returns to a near-circular orbit at the end of its first Kozai cycle. The 
system periodically crosses the secular separatrix whenever e2 goes below or 
above ~ 0.1. Once the system is trapped in the tight apsidal libration regime, 
the eccentricity ratio is maintained at ~ 10, thus ei also oscillates together with 
e2 with a smaller amplitude. 



3. Discussion 

The rigid rotation of the planetary orbits and synchronized eccentricity oscil- 
lations are frequently observed among our simulations. The criterion for the 
up-libration in Equation [3] can be easily satisfied if the Kozai oscillation ampli- 
tude of the outer planet is sufficiently large. However, the system also needs to 



3 Here we assume that the Kozai oscillation of the inner planet is suppressed by the secular 
torque applied by the outer planet. Normally this assumption is valid for a large range of a± 
when mi < m.2. 
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maintain its stability for sustain able Kozai oscil lations. One necessary stability 
requirement is the Hill stability (jGladman 19931 ). The Hill stability needs to be 
maintained at all times including when e2 is near its maximum. Another case 
of instability arises when the initial conditions are chosen such that the sys- 
tem undergoes down-libration (Eq. [2|) is satisfied. During down-librations the 
apsidal vectors are anti-aligned (Aw ~ 180°), and e2/e\ ~ 1 or smaller. If the 
down-libration is maintained until the outer planetary orbit embarks on the first 
Kozai cycle, then e\ is excited to a value larger than e2- Such an over-excitation 
leads to a fairly chaotic behavior of the inner orbital eccentricity and will most 
likely result in collision of the inner planet with the central star or ejection from 
the system. 

The synchronized eccentricity oscillation mechanism offers a few theoretical 
implications. First, in this mechanism the outer planet helps the angular mo- 
mentum transfer from the companion star to the inner planet. Propagation of 
eccentricity disturbance by a stellar passerby through a system of planets has 
been previously studied by IZakamska Tremaine ( 20041 ) . In the case of two 
planets in a binary, the cyclic propagation of the stellar perturbation through 
the outer planet can effectively excite the inner planet's eccentricity to oscil- 
late. This eccentricity propagation mechanism may imply the presence of an 
undetected secondary planet around known single-planet systems in binaries. 
Especially for single planets in wide binaries with the Kozai period too large 
compared to the estimated age of the system, yet with non-negligible observed 
eccentricities that are difficult to be explained by other dynamical processes, 
there is a possibility that there is an undetected secondary planet propagating 
the eccentricity excitation from the distant companion. 

Another dynamically relevant outcome is the subsequent evolution of the 
planets after the Kozai mechanism induces instability to the system . Many 
systems in our simulations have suffered close encounters between planets that 
resulted in losing a planet from the system by ejection or collision. The or- 
bital radius and the eccentricity of the survived planet can take wide ranges 
of va l ues, as investigated by planet s cattering experiments (jJuric Tremaine 
20071 : IChatteriee. Ford fc Rasiol l2007fl . The survived planet may resume the 



Kozai mechanism, but now the planet's orbital eccentricity is likely to have 
been already excited to a significantly larger value after ejecting the other 
planet. Since the fraction of the time a planetary orbit remains highly eccen- 
tric during the Kozai oscillations is strongly dependent on its initial eccentricity 
( Holman. Touma Sz Tremaind 19971 ). single-planet systems in binaries with an 
episode of early scattering has a much higher probability to be observed at large 
eccentricities. For example, known single planets in binaries with very large 
orbital eccentricities (e.g., HD 80606 b, e = 0.93 and HD 20782 b, e = 0.92) 
strongly suggest the past occurrence of planet ejections in these systems, en- 
abling the eccentricities of the survived planets to oscillate at systematically 
larger values. 
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